Decoherence Correction by the Zeno Effect and Non-Holonomic Control 
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We show that multidimensional Zeno effect combined with non-holonomic control allows to effi- 
ciently protect quantum systems from decoherence by a method similar to classical coding. Contrary 
to the conventional approach, our method is applicable to arbitrary error-inducing Hamiltonians and 
general quantum systems. We also propose algorithms of finding encoding that approaches the Ham- 
ming upper bound along with methods of practical realizations of the encodings. Two new codes 
protecting 2 information qubits out of 7 and 4 information qubits out of 9 against a single error 
with arbitrarily small probability of failure are constructed as an example. 

PACS numbers: 03.65.-w, 03.65.Fd, 03.67.Lx, 32.80.Qk 



Decay and decoherence are the two processes that oc- 
cur in open quantum systems. Two difïerent relaxation 
times T\ and characterize these processes in two-level 
systems. Though both processes represent irreversible re- 
laxation, they are of a rather different nature. The num- 
ber of occupied states given by the rank of the density 
matrix - an analog of the phase volume in the classical 
case - may change in an elementary act of relaxation, 
and usually does it towards diminishing. By contrast, 
this number does not change in an elementary act of de- 
coherence, but yet increases after being averaged over 
subsequent elementary acts. 

In this Letter, we show how one can completely sup- 
press decoherence in a universal way. We employ bàsic 
ideas of classical error-correcting codes and their powerful 
extensions [1] that make use of extraspecial and Clifford 
groups to build quantum codes from their classical coun- 
terparts. However, our approach considerably deviates 
from this method, since it relies on the full group of uni- 
tary transformations in the entire Hilbert space. We also 
employ the idea known as the Zeno effect [2] that states 
that periòdic measurements allow to keep a quantum sys- 
tem in the initial state for as long as needed, provided 
that the period T is short relative to the relaxation time 
T2. However, contrary to the conventional Zeno effect, 
we do not intervene into the entire quantum system. In- 
stead, we only affect the auxiliary part, the ancilla, to 
protect the quantum state of the main part. As a gen- 
eral tool for protecting any system from decoherence, we 
suggest the non-holonomic control [3], which allows one 
to perform any predetermined transformation including 
encoding in the Hilbert space of a quantum system. 

An important example emerges when a classical binary 
code is compared with its conventional quantum counter- 
part [1,4]. Both codes add the auxiliary part to increase 
minimal discrepancy between possible states. However, 
classical codes use a finite alphabet mapping k informa- 
tion bits onto a single code vector of a bigger length n. By 



contrast, the state of k qubits is a continuous complex- 
valued vector in the entire 2 fc -dimensional Hilbert space 
Hk ■ When the original system is expanded by n — k an- 
cilla qubits in some state |5), the extended state vectors 
form a subspace T~í k ~ — |S) (a\Tí n in the bigger space 

Tí n . Quantum encoding C is a unitary transformation 
that maps the subspace Ti. k ~ onto a 2 fc dimensional code 

subspace C — CH k ~ in TL n . Note that C spans over all 

columns of the rectangular matrix C\a), and has zero 
projections on all 2" — 2 k vectors that form the columns 
of any other matrix C \a) with \a) _L \a) in the ancilla. 

There is also a difference in the decoding process. A 
corrupted codeword can be fully retrieved and corrected 
by a classical code. By contrast, in quantum mechanics, 
retrieval of the entire state vector by a measurement de- 
stroys the data. A natural way of quantum correction is 
therefore to find an encoding C and decoding C -1 that 
can move all possible errors to the ancilla. In this case, 
resetting the ancilla to the state \a) gives a projection of 
the entire state vector exactly to its original position in 
the subspace 7í k ~. To this end, given any set of errors 

j-E m j, we should find a code C ={\v)} such that all pos- 
sible error vectors Ei \v) are orthogonal to C and to each 
other [5] according to the condition 
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Formally, this requirement coincides with the definition 
of a code C (see [4] p. 436) written for zero trace matrices 
E m [6] , although the errors we consider are not restricted 
to the extraspecial group in contrast to the original for- 
mulation. It can be seen that Eq.(l) also corresponds to 
the classical Gilbert-Varshamov [7] bound. 

Below we describe a recipe for constructing codes that 
are "nearly" orthogonal to error vectors with any prede- 
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termined accuracy, so that the corrupted state vector falls 
arbitrarily close to the original code vector, after being 
projected on the code subspace. It allows one to replace 
the rather strict condition Eq.(l) by a weaker condition, 



of the ancilla is given by the commutator 



(v'\È m \v) =0; V \v) , \v') e C; ~ÍE m e {È m } 
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Note that the latter can be thought of as an analog of the 
classical Hamming bound [7] . We emphasize that our ap- 
proach extends beyond the quantum systems composed 
of identical two-level partides and is vàlid for any quan- 
tum system separable into a main part and an ancilla. 
In fact, this approach only depends on the total numbcr 
M of the crror-inducing Hamiltonians |s m | regardless 
of their specific matrix structure. 

To design codes that meet the conditions of Eq.(2), we 
will use the Zeno effect. The regular Zeno effect allows 
one to restore the initial state vector up to the second- 
order terms in T/T2, which are very small provided that 
T <C T 2 . The essence of this phenomenon is based on 
the heurístic that within a short time interval, a unitary 
quantum evolution moves any state vector in a direc- 
tion orthogonal to the vector itself. In our version of the 
Zeno effect, the subsequent projection is done to the 2 k - 
dimensional subspace 7~L k ~ instead of the state vector. 
Thercforc, our primary goal is to find a code C such that 
for short timc intervals, any error Hamiltonian in |-E m | 
moves any code vector in a direction perpendicular to C. 

To implement the Zeno effect, consider the state vec- 
tor \S) = \s) <g> \a) of the compound system formed by 
the information part in the state \s) and the ancilla in 
the state \à). The vector \S) undergoes an uncontroUed 
unitary evolution 



M 
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iE m Jf m (t)dt 



M 



T-iJ2 E m fm(t)dt 



2=1 



(3) 



where the time ordering of the product is implícit. To il- 
lustrate the idea of the approach, we simply treat the 
errors E m as Hamiltonians of interactions with exter- 
nal random fields f m (t) that produce an uncontroUed 
evolution of the system. We also assume that over 
the Zeno period T the diffcrcnt fields have actions 



i(x)dx 



<C 1 that are so small that only the 

identity operator / and the first order terms ~ Jf m (t)dt 
are important in the Taylor series, whereas the highcr 
orders can be ignored. 

Let p sc = \s) <S> (s| be the density matrix of the main 
part before the action of errors. Then p = \S) <& (S\ — 
\a) p sc (a\ is the density matrix of the entire system. The 
variation Sp sc = {à\5p\a) of the density matrix of the 
main system after the perturbation Eq.(3) and resetting 
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Therefore p sc satisfies the master equation 



.dp s 



M 



dt 



h e , p sc ; h e =^2f m (a\C 1 E m C\a) 
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with effective Hamiltonian h e . Then requirement Eq.(2) 
implies that h e = and p sc =const, whereas the relation 
\v) = C \s) <g> \a) holds for the code vectors and the states 
of the system. 

Our next step is to design codes that satisfy Eq.(2). 
Let N — 2 k and A = 2 n ~ k be the dimensions of the 
Hilbcrt spaces formed by the main part (information sub- 
system) and ancilla, respectively. Then Eq. (2) is a set 
of M x TV 2 equations taken over A x iV 2 matrix elements 
(a| (g> (s'\C\s) <S> |5) of the encoding operator that de- 
fine the rectangular TV x NA matrix C\à). To solve this 
system, we wish to limit the number of equations by the 
number of free variables. This gives the Hamming upper 
bound in the form M < A. To specify Eq. (2) further, 
let the states \s) (g> |5) correspond to the first N positions 
of the state vector, which implies that each of the M 
matrices C~ 1 E m C has all-zero upper left N x N còrner. 
Standard methods of linear àlgebra do not give a recipe 
of finding a linear transformation C that simultaneously 
sets to zero the còrners of M different matrices. The first 
result of the paper is an iterative algorithm which docs 
so, and hence finds a code C that satisfies the conditions 
of Eq.(2) given M arbitrary error Hamiltonians E m . 

To illustrate the method, we first find a single vector 
\x) of length N x A that is orthogonal to all M error 
vectors E m \x). In other words, \x) should satisfy M 
quadratic equations 



(x\ E m \x) = 0. 



(6) 



Our input is a randomly selected vector \x) of unit length. 
If, by chance, \x) is orthogonal to all vectors E m x), the 
problem is solved. Otherwise, there exists a linear com- 
bination X)m=i lmE m \x) that minimizes the length of 
the vector X)m=i lmE m \x) + \x). To find it, we employ 
the Standard variational principle and solve the corre- 
sponding set of linear equations for j m . In our further 
iterations, we set \x) — ► J2m=i \lm,E m \x) + \x), then 
normalize the new vector to unity, and repeat the pro- 
cedure until it converges. The same algorithm serves to 

find the entire encoding matrix C. For this purpose, one 

— — > 

should simply consider a supervector x — C \ a) of length 
N 2 A by appending the columns of C\a). Then the or- 
thogonality conditions (vi \v s ) = taken for the different 
columns s,l, s ^ l, and the conditions {vi\ E m \ v s ) = 
give independent equations for this supervector. 
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This procedure has shown a ràpid rate of convergence 
for all the examples considered, including two new singlc- 
crror correcting codes with (n, k) parameters (7, 2) and 
(9,4) [8]. Note that these codes are designed in the 
Hilbert space and pcrform singlc error correction with 
any predetermined accuracy that depends only on T. 
However, neither of the two codes can exist for conven- 
tional quantum techniques that satisfy Eq.(l). More gen- 
erally, our conjccture is that the above procedure con- 
verges if the length of a supervector x is no less than the 
mrmbcr of quadratic equations Eq.(6), or, equivalently, if 
the Hamming bound M < A holds. For this conjecture, 
note also that all matrices E m in Eq.(6) have zero traces. 

Once a code C is found, a question arises as to how the 
system can be moved into this subspace. Thus, we also 
nccd to dcsign an cfhcient encoding C. For the codes 
built on the Clifford group, this can be done by ap- 
plying a polynomial number of Standard operations like 
Hadamard gates. It is not the case for an arbitrary quan- 
tum system, since for genèric error sets |_E m |, C is also a 

genèric unitary matrix, and hence an alternative mcthod 
of coding is required. One can achieve this feasibility 
with the help of non-holonomic control, by applying a 
sequence of two different natural Hamiltonians H\ and 
Hi that satisfy certain conditions [3] . The second result 
of the paper is an algorithm that allows one to determine 
a set of non-holonomic control timings t\, t 2 , t 3 , ■ ■ ■ ,tM' 
of a reasonably short length M' resulting in the encoding 
transformation 

(J — e ~itM' H 2 . . . e -it 3 H le -it 2 H2 e -it 1 H 1 ^ 

The number M' equals or slightly exceeds the total num- 
ber MN 2 of rclevant error matrix elements. 

The key idea of the method is to incorporate the search 
for the non-holonomic control parameters ti into the it- 
erative algorithm for the constants j m . Being written in 
the form of an N x N matrix j m , these constants relate 

the supervector C\a) to its change Em=i^mC|5)7 m , 
thus allowing us to find the first-order variations 5 m — 
(5| C^Í^C |5) 7 m +j m (a\C- 1 È m C\a) of the error 
matrices projected on the code space. We therefore 
choose the variations Sti of the control timings such that 
the corrcsponding increments of the error matrices are 
equal to 6 m for each to, that is, 

(à\ C- x E m C \a) 7m + jl (à\ C-'ÈlC \a) 




Then we continue the iterations, while varying the big 
parameter (3 to gain maximum convergeney. The num- 
ber M' of the parameters ti can be chosen slightly above 
the minimum number MN 2 in order to avoid a prema- 
ture stop if the MN 2 x MN 2 supermatrix in the paren- 



theses of Eq.(8) becomes singular for some intermediate 
approximation C of the coding operator. 

This error-correction method becomes asymptotically 
exact for T — > 0. However promising for small quan- 
tum systems, it has complexity that is exponent in the 
number of partides n. Therefore, we also propose an 
alternative strategy applicable to large systems, which 
is the third result of the paper. If the computational 
time required to find the MN 2 x MN 2 matrix becomes 
prohibitive, thus hindering any way to exactly meet the 
conditions of Eq.(2), one can minimizc the rate of er- 
ror aceumulation in the density matrix, by minimizing 
the effective Hamiltonians in Eq.(5). For the systems 
with binary interaction an exponential reduction of the 
error-aecumulation rate can be achieved if the projection 
rate 1/T grows faster than n 2 . Moreover, the code rate 
R = k/n tends to unity as 1 — o(n 2 T). 

In particular, one can reduce all effective Hamiltonians 
h m = fm (5| C~ l E. m C \à) by applying a random encod- 
ing C 7 which allows to inhibit the rate of decoherence 
by the order of ancilla dimension A = 2 n ~ k . Indeed, 
for individual qubits, a typical single-particle or binary 
interaction E m is given by a sparse matrix that strongly 
couples each state to sufficiently few other states. By con- 
trast, most matrix elements are zeros. A genèric unitary 
encoding C smears out this coupling over all 2" states of 
the entire system, thus making a typical matrix element 
of C~ 1 E m C 2™ times smaller than the typical nonzero 
matrix element of E m . Since only 2 k of these states 
contribute to the projected interactions h mi the overall 
gain is equal to 2"~ fc . From a more general viewpoint, 
this random encoding rests on the fact that the decoher- 
ence rate is a polynomial in the number n of partides, 
whereas the dimension of the entire Hilbert space is an 
exponent in n. It allows to store quantum information in 
strongly entangled states of many-body systems that are 
only weakly affectcd by physically realistic perturbations. 

This random encoding also allows to specify a relation 
between the dimension N — 2 k of the information sys- 
tem and the maximum number t of possible errors caused 
by binary interactions. Given K single-particle quantum 

states, there exist M — J2l=i (?) {K 2 ~ l) different pos- 
sible error matrices (here K — 2 for qubits). Since the 
energy of binary interaction seales as n 2 , each error ma- 
trix yields a decoherence rate of W — n 2 /T 2 . This rate 
is reduced A — 2 n ~ k times after the matrix elements are 
smeared out over the entire system. Thus, given n — > oo 
identical partides, one can correct t errors with high fi- 
delity if the total decoherence rate W2 k ~ n M vanishes. 
This condition W2 k -' n M -> again gives the Hamming 
bound [7] 

1 - fc/n-log 2 M 1/n > 0. (9) 

Hence, to asymptotically meet this classical bound, we 
can randomly choose an (n, fc)-code thus leaving the error 
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patterns of weight t or lcss uncorrected with a vanishing 
probability. 

Finally, note that the sequential application of two 
different Hamiltonians suggested by Eq.(7) rcsults in a 
completely genèric transformation, provided that: (i) the 
numbcr M' of the control timings í; is of the order of 
log 2 (N x A); (ii) the operators l·L\ and Hi satisfy the 
conditions of non-holonomic control [3] ; and (iii) the time 
intervals í; are long enough to ensure big acquired actions 
11-^1,2 3> fi. Thcrcfore, given a system, the main prob- 
lem is to find two controllable interactions that satisfy the 
conditions needed for non-holonomic control, and allow 
to perform the decoding transformation C^ 1 at the same 
level of complexity. 

We present such a pair H\ , H 2 for the example of the 
molecule proposed in [9] as a toy quantum computer of 
7-qubits, for which the applicability of both the exact 
and the approximative approaches was explicitly verified 
by determining the control timings and the reduction of 
cffcctivc coupling, respectively. The first operator is the 
nuclear magneto-dipole interaction 

H 1= J2B xl·lj a^ (10) 

3 

of the molecule in the ground ro-vibronic state \g) with 
a static external magnètic ficld B x oriented along the 
x-axis. The second operator is the Raman coupling 

H2 = Y1 E ^-vgw^®^ 1 

i<j a,b=x,y,z (11) 

in a field B R oscillating at a high-frequency u> detuned 
by Suj — Lütr — ui from the frequency uitr of transition to 
an exited ro-vibronic state |e). Here fj,j = (g\ Hj(r) \g) is 
the average value of the coordinatc dependent gyromag- 
netic ratio fij (r) of the j-th nucleus in the ground state, 
and fi[ 2 j = (g\ ^i(r) |e) (e| fJ-j(r) \g) represcnts the matrix 
element of the Raman transition. By changing the di- 
rection of B x from positive to negative and reversing the 
detuning Sui, one can vary the signs of the Hamiltonians, 
which allows one to easily construct the inverse transfor- 
mation C~ x just by inverting the order of the sequence 
ÍM', ■ ■ ■ ,ti hi Eq.(7) along with the change of the signs. 
We do not present here [8] a cumbersome set of M' <~ 

2 5+2x2 

non-holonomic control timings found explicitly 
with the help of the approximative algorithm for the 
above single error-correcting code that protects two out 
of seven qubits in the toy computer. Nor do we give the 
explicit matrix elements of the effective error Hamiltoni- 
ans Eq.(5) obtained as a result of ten sequential appli- 
cations of H\ and Hi that reduce the decoherence rate 
of the two selected qubits. We only note that our nu- 
mèric results are completely consistent with the given 
estimates, whereas the absolute vàlues of the matrix ele- 
ments obey a Maxwellian-like statistics. 

Summarizing the main rcsults of the paper, we con- 
sider open quantum systems that consist of two parts - 



the main part and an ancilla. The latter is being peri- 
odically reset to a given quantum state while the main 
part is protected against uncontrolled decoherence with 
the aceuracy proportional to the period. The numbcr 
of different error Hamiltonians should not exceed the di- 
mension of the ancilla Hilbert space. Initially disentan- 
gled, the two parts are driven to a highly cntanglcd state 
by means of the unitary non-holonomic control. This 
new state belongs to a code space such that any error 
can only move the code vectors in a direction perpendic- 
ular to this space. The inverse unitary transformation 
brings the system to a state that differs from the initial 
state only in the ancilla part. After resetting the an- 
cilla, the system returns back to the original state. We 
propose new algorithms that give the corresponding code 
spaces and parameters required for non-holonomic con- 
trol. Thcsc algorithms have complexity proportional to 
the dimension N of the entire Hilbert space. For large 
systems, the computational time required to find a code 
space becomes prohibitive. However, a coherent control 
sequence needs only logTV switches to strongly reduce 
the decoherence rate and bring the system to a highly 
entangled state weakly sensitive to the decoherence. 
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